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We determine by small-angle X-ray scattering the structure factor of hydrophobic particles inserted within 
lamellar surfactant phases for various particle concentrations. The data is then analyzed by numerically solv¬ 
ing the Ornstein-Zernicke equation, taking into account both the intra- and inter-layer interactions. We find 
that particles within the same layer repel each other and that the interaction potential (taken as independent 
of the concentration) has a contact value of 2.2 ksT and a range of about 10 A. If the amplitude is allowed 
to decrease with increasing concentration, the contact value in the dilute limit is about 5 ksT, for a similar 
range. 
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I. INTRODUCTION 

The elucidation of membrane-mediated interaction be¬ 
tween inclusions in the cell membrane (such as integral 
proteins or membrane-active antibiotic peptides) is of 
paramount importance for understanding their biological 
activity. Since the underlying problem is the organization 
of the membrane -seen as a two-dimensional complex 
system- it has become clear over the last three decades 
that the concepts developed in soft matter physics for the 
understanding of self-assembled systems are operative in 
this context and that ‘simple’ models can yield valuable 
information. 

It is therefore not surprising that sustained theoret¬ 
ical efforts attempted to provide a detailed description 
of these complex systems; they are either continuum- 
elasticity theoriesi^— or more detailed models taking into 
account the molecular structure of the lipid bilayer— 

In the few cases when the models were validated and 
refined using experimental data, these approaches were 
often successful. For instance, the lifetime of the gram¬ 
icidin channel is known to depend on the thickness^ 
and tensioniS of the membrane; it was shown that a 
continuum-elasticity model provides a satisfactory de¬ 
scription of this phenomenoiiii. 

However, this body of theoretical work has not yet 
been matched by the experimental results, the first of 
which were obtained by directly measuring the radial dis¬ 
tribution function of membrane inclusions using freeze- 
fracture electron microscopy (FFEM)i2i“— . These data 
were compared to liquid state models and could be de¬ 
scribed by a hard-core model with, in some cases, an 
additional repulsive or attractive interactioni^^— . How¬ 
ever, FFEM has not been extensively employed since, cer¬ 
tainly due to the inherent experimental difficulties; fur¬ 
thermore, the particle distribution observed in the frozen 
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sample is not necessarily identical to that at thermal 
equilibrium. 

A step forward was taken by Huang and collabora¬ 
tors, who showed that the distribution of membrane in¬ 
clusions within the plane of the layers can be studied us¬ 
ing small-angle X-ray and neutron scatterin g^^i^° . First 
of all, these techniques are perfectly adapted to the typi¬ 
cal length scales to be probed. The measurement is aver¬ 
aged over a large number of particles and over long times, 
without perturbing the system; thus, one has access to 
the structure factor of the interacting particles. How¬ 
ever, in these studies, at most two peptide-to-lipid con¬ 
centrations PjL were investigated for each system and no 
values were given for the interaction potential, the data 
being explained in terms of purely hard-core interactions. 

Building upon this work, we recently studied systems 
with a varying density of inclusions. This is indispens¬ 
able since, even though each structure factor (taken sep¬ 
arately) can be described by a hard-core model, the 
apparent radius obtained changes with the concentra¬ 
tion, signalling the presence of an additional interac¬ 
tion. For instance, we were able to measure the in¬ 
teraction potential of alamethicin pores in dimyristoyl- 
phosphatidylcholine (DMPC) bilayers. We showed that, 
aside from the expected hard-disk repulsion (with a ra¬ 
dius corresponding to the geometrical radius of the pore) 
the pore interaction exhibits a repulsive contribution, 
with a range of about 3 nm and a contact value of 
2.4 /cbT—. On the other hand, for gramicidin pores in 
dilauroyl-phosphatidylcholine (DLPC) bilayers, while the 
interaction is still repulsive, the parameter values are 
quite different, with a higher contact value and much 
shorter range. In the latter case, the interaction was 
shown to decrease with the pore concentration, in agree¬ 
ment with the hydrophobic matching model^^. These 
results are in qualitative agreement with recent theoret¬ 
ical predictions^i^^ and can be used as a test for other 
theoretical and numerical results. 

The biological molecules cited above act specifically 
by insertion within cell membranes, so their affinity for 
lipid bilayers is assured, which is a significant advantage. 
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However, their use as membrane probes has important 
shortcomings: their X-ray scattering contrast is low and 
the positions of the constitutive atoms are not always 
well-defined (membrane proteins can adopt various con¬ 
formations as a function of the environment, while the 
pores formed by antimicrobial peptides often comprise 
variable numbers of monomers). These features con¬ 
tribute to the difficulty of obtaining high-quality scater- 
ing data on membranes with inclusions. 

Fortunately, if one is interested in the general physical 
properties of membranes seen as two-dimensional com¬ 
plex fluids, rather than in the behaviour of a specific 
active molecule, other inclusions can be used. For in¬ 
stance, we have shown recently that tin oxo-clusters can 
be inserted within surfactant bilayers and that their in¬ 
teraction potential can be determined^^. 

These types of hybrid organic-inorganic particles have 
the following advantages: 

• They are perfectly monodisperse (actually iso- 
molecular) and “rigid” (their atomic configuration 
is well-defined). 

• Their scattering contrast is high (due to the pres¬ 
ence of metal atoms). 

• Their surface properties can be tailored by chang¬ 
ing the nature of the peripheral groups. 

The purpose of this article is to develop a full anal ysis 
of the system, beyond the initial approach used in Ref. M. 
Two main improvements can be noted: 

• We were able to determine the complete structure 
factor, S{qr,qz), while the initial study was only 
concerned with the equatorial slice, S{qr,qz = 0). 
This required new experimental data, for q^ ^ 0. 

• We employed a more elaborate statistical model 
(based on the numerical solution of the Ornstein- 
Zernicke equation with the Percus-Yevick closure) 
to relate the structure factor to the interaction 
potential; our first approach was based on the 
random-phase approximation (RPA). 

We have also checked the reproducibility of the results 
(see Appendix lAT). by comparing the old and new data, 
thus verifying both the reliability of the method and the 
stability of the samples. 


II. MATERIALS AND METHODS 


The membranes were composed of 
dimethyldodecylamine-N-oxide (DDAO), a single¬ 
chain zwitterionic surfactant. The DDAO (purchased 
from Sigma-Aldrich) was first dried in vacuum for 20 h 
to remove any residual water— and then dissolved in 
isopropanol. 

The two stock solutions were then mixed to yield the 
desired BuSnl2/DDAO ratio and the mixtures dried in 
vacuum, yielding a final mass of about 200 mg for each 
sample. Varying amounts of water were then added so 
that the mixtures were in the fluid lamellar La phase^^. 
For DDAO, the molecular weight is 229.40, the density 
is 0.84 g/cm^ and the bilayer thickness is 25 ± 1 
resulting in an area per surfactant molecule Addao = 
37.8 A^. The molecular weight of the BuSnl2 is 2866.7, 
with a density of 1.93 g/cm"^—. The (two-dimensional) 
number density of particles in the plane of the membrane, 
n, is calculated using the data above and neglecting the 
possible increase in bilayer surface due to the inclusions. 

The lamellar phases were then drawn into flat glass 
capillaries (VitroCom Inc., Mt. Lks, NJ), 100 fim thick 
and 2 mm wide by aspiration with a syringe and the cap¬ 
illaries were flame-sealed. Good homeotropic alignment 
(lamellae parallel to the flat faces of the capillary) was ob¬ 
tained by thermal treatment, using a Mettler FP52 heat¬ 
ing stage. The samples were heated up to the isotropic 
phase (at 130 °C) and then cooled down to the lamellar 
phase at a rate of 1 °C/min. 

We studied the samples by small-angle X-ray scatter¬ 
ing (SAXS), on the the bending magnet beamline BM02 
(D2AM)^ of the European Synchrotron Radiation Facil¬ 
ity (ESRF, Grenoble, France). The photon energy was 
11 keV and the sample-to- detector distance about 27 cm, 
with a scattering vector in the range 0.04 < q < 0.9 A 
The detector is a charge-coupled device Peltier-cooled 
camera (SGX90-1300, from Princeton Instruments Inc., 
NJ) with a detector size of 1340 x 1300 pixels. Data pre¬ 
processing (dark current subtraction, flat field correction, 
radial regrouping and normalization) was performed us¬ 
ing the bm2img software developed at the beamline. The 
scattering geometry is discussed in section HVl 


III. MODEL AND ANALYSIS 

Since the particles are identical in shape, the scattering 
intensity can be written as the product of a form factor 
(only depending on the internal constitution of the par¬ 
ticles) and a structure factor, which describes the inter¬ 
action between particles^: I{q) = S{q) ■ \F{q)\^, with: 


We used as inclusions hybrid nanoparticles consist¬ 
ing of a tin oxide core decorated with butyl chains. 
The complete formula is: {(BuSn)i20i4(0H)6}^+(4- 
GH3C6H4S0;^)2, shortened to BuSnl2 in this pa per The 
synthesis and structural details are given in Ref. The 
BuSnl2 particles were dissolved in ethanol prior to use. 


S{qz,qr) 


1 
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^ exp {-iqvk) 


/c=l 


( 1 ) 


where N is the number of objects and object “0” is taken 
as the origin of the coordinates. 
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We expect the form factor of the particles to be dom¬ 
inated by their inorganic core, as the electron density of 
the butyl chains is similar to that of the dodecyl chains 
within the bilayers and to that of ethanol. This assump¬ 
tion is confirmed by the intensity at higher scattering vec¬ 
tors (g > 0.5 A ), which is well described for all samples 
by the form factor of a sphere \Ff{R, g)p, with a radius 
(used as a free fitting parameter) R = 4.5 ± 0.2 A, in 
good agreement with the average radius of the tin oxide 
“cage” estimated from the crystallographic data. 

Dividing the measured intensity by the form factor 
above yields the structure factor S{q). Standard error 
propagation then yields the uncertainty, but the resulting 
values (corresponding to the statistical error) are much 
smaller than the discrepancy between the fitting model 
and the experimental data (the goodness-of-fit function 
per data point is much larger than one). Recognizing 
that this discrepancy is mostly due to systematic errors 
in data acquisition, background subtraction etc. and also 
possibly to the inadequacy of the theoretical models, we 
assign to each data point a fixed uncertainty cr; in the 
following, we take for definiteness a = 0.05, which is a 
rough estimate of the difference between model and data, 
insofar it yields values of unit order. 

Determining the structure factor for systems composed 
of many particles interacting via a known potential is one 
of the fundamental endeavors of liquid state theory. 

Two cases must be considered: 

• The situation when the particles are dissolved in 
a simple solvent will act as a “reference state”; 
more precisely, it will allow us to check whether 
the “naked” particles can be described simply by a 
hard core repulsion or whether an additional term 
must be considered and what its characteristics are. 

• The most interesting situation is of course that 
of particles inserted within the surfactant bilayers, 
which will be the main result of our study. Due to 
the anisotropic nature of the matrix, the analysis 
must account for the difference between the inter¬ 
action between particles contained within the same 
layer and that from one layer to the next. 


A. Structure factor of an isotropic solution 

In solution, the interaction between particles is taken 
as isotropic, described by a two-body potential depend¬ 
ing only on the distance between their centers: V 3 d(t’) = 
Vsnir)- Consequently, the structure factor S{q) is 
equally isotropic, depending only on the magnitude of 
the wave vector q. 

Using as input the interaction potential, the structure 
factor was determined numerically by solving iteratively 
the Ornstein-Zernicke equation with the Percus-Yevick 
closure condition. The algorithm is based on the method 
of Lado ^^i^^ ; see Appendix |B] for the details. The in¬ 
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FIG. 1. Top: Experimental three-dimensional structure fac¬ 
tors for BuSnl2 particles in ethanol (symbols) and fits with 
the Gaussian model discussed in the text (solid lines) for five 
volume concentrations (j>. S{q) decreases with increasing 0. 
Bottom: Interaction potential obtained from the data above 
via the numerical procedure discussed in Section IIII Al (solid 
line and shaded area for the uncertainty range) and within 
the RPA approximation (solid line and symbols). For com¬ 
parison, we also show the RPA results obtained in Ref. 
(dashed line). 


teraction potential is described by a hard core with di¬ 
ameter 2R = 9.4 A and an additional “soft” component, 

modelled either as a Gaussian centered at the original: 
2l 


V{r) = f/oexp 


t(i)' 


r > 2R or as a decreasing 


exponential with contact value Uc and decay range 
U(r) = Ucexp [—(r — 2i?)/^] r > 2R. 
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B. Structure factor in the lamellar phase 

The phase (and hence the structure factor) is now 
anisotropic. We assume that there is no in-plane order¬ 
ing (the particles form a 2D liquid in the plane of the 
bilayer), so that S only depends on the absolute value of 

the in-plane scattering vector qi + Qy and on the 

scattering vector along the membrane normal q^. 

For an ideal gas (no interaction) S = 1, while if the 
particles interact only in the plane of the bilayer and not 
from one bilayer to the next, S = S{qr) only. In the 
general case, the structure factor can be developed as^S: 


S{qr,qz) = So{qr) + 2 ^ cos{qzdm)Smiqr) ( 2 ) 

m—1 

with each partial structure factor Sm describing the in¬ 
teraction between particles situated m layers away: 

pCO 

Sm{qr) = Som + 2np rdr Jq( rg^) [ffm(»") “ 1] (3) 

Jo 

where m > 0 and gmif) is the (normalized) probability 
of finding a particle at a distance r in bilayer m, knowing 
that a particle is present at the origin in bilayer 0. 

This expansion merely reflects the discrete nature of 
the stack along the z direction. It is especially conve¬ 
nient since generally the interaction does not extend very 
far along z. Note also that Equation ([2]) is only valid for 
the geometrically perfect case when the distance between 
neighboring layers is rigorously d. In practice, this dis¬ 
tance varies due to thermal fluctuations and to frozen-in 
defects, leading to a smearing of the diffraction pattern 
at high g^ values. We describe this effect phenomenolog¬ 
ically by a Lorentzian factor. Finally, we use: 


S{qr,qz) = Soiqr) + (g,) (4) 

1 -b (qzCJr 

where the disorder parameter a has units of length. For 
thermal fluctuations, tr^ = ((zi — zq)'^) cz 77(d/7r)^, where 
77 is the Caille parameter. 

As we will see below, the structure factor Si{qr) due 
to nearest-bilayer interaction has a localized peak. For 
convenience, we describe it by a Gaussian function: 


Si = Ai exp 


{q ^max) 
2Ag2 


( 5 ) 


The partial structure factors So{q) and 5'i(g) are simi¬ 
lar to those describing a two-dimensional binary mixture 
AB, with the formal identification So = Saa = Sbb and 
Si = Sab- As we will see below, these two components 
are enough to describe the experimental data, so there is 
no need to go beyond m = 1 (interaction between adja¬ 
cent bilayers). 


The (experimentally determined) functions So{q) and 
Si (g) can be described -using well-established results in 
liquid state theory- in terms of the interaction potentials 
Vo('r) (between particles within the same layer) and Vi(r) 
(between adjacent layers). The derivation of the model 
and the implementation details are given in Appendix IbI 


IV. RESULTS 
A. Interaction in solution 

The 3D structure factors of BuSnl2 particles in ethanol 
are shown in Figure [T] (top) for five different concentra¬ 
tions (lines and symbols). They are fitted with the Gaus¬ 
sian interaction model discussed in Section ITlI Al yielding 
Uq = IBksT and ^ = 5.1 A, such that the potential at 
contact is Uc = 3.5 kBT. The fits are shown as solid 
lines. The corresponding interaction potential is shown 
in Figure[T] (bottom) as solid line. The shaded area repre¬ 
sents the uncertainty (see Appendix IB 31 for details). For 
comparison, we show in the same figure the potential ob¬ 
tained in the RPA approximation for the current data 
(solid line and symbols) and for that in Ref. [13 (dashed 
line). The exponential model yields very similar results, 
in terms of goodness of fit and in potential shape for 
parameter values Uc = 4.5 fcsT and ^ = 2.1 A. 


B. Interaction between layers 

In order to quantify the interaction between the lay¬ 
ers (along the smectic director z), one needs access to 
the complete structure factor S{qr,qz)- This is achieved 
by using the experimental configuration described in Fig¬ 
ure m as first discussed in Ref. the incoming X-ray 
beam (with wave vector ki) is incident upon the flat cap¬ 
illary at an angle a to its normal n, which coincides 
with the smectic director z. A point on the 2D detec¬ 
tor uniquely defines an outgoing wave vector kg for the 
scattered signal, and thus a scattering vector given by: 
q = kg — ki. One can therefore assign to each pixel values 
for g^ = g • z and qr = \q- qzz\- 

The raw scattering image is then re-gridded to an ap¬ 
propriate region of the reciprocal space, {qr,qz)- For each 
pixel in the target space, the algorithm identifies the cor¬ 
responding point in the starting image, and the closest 
9 pixels are averaged to yield the final intensity value. 
It is easily shown that the accessible range in reciprocal 
space is triangle-shaped, with an angle 2 q: at the origin. 
Finally, the resulting image is divided by the form factor 
of the particle to yield the structure factor. 

As an illustration, we show in Figure[2]re-gridded data 
at four different incidence angles a for the sample with 
n = 1.97610“^ A . A modulation along g^ is clearly 
visible, indicating the presence of the Si component and 
hence of an interaction between layers. The images also 
display a linear slope in g^, probably due to imperfect 
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FIG. 2. Left: Sketch of the experimental conhguration. Right: Re-gridded scattering data for various values of the incident 
angle a (indicated within the images). 
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FIG. 3. Gomparison between the experimental data for 
S{qr,qz) — S{qr,qz = 0) derived from the image taken at 
q: = 60° in Figure [2] and the model ([4]). 


background subtraction. This artefact is removed prior 
to further treatment. 

Although the physically relevant functions are S'o(9r) 
and Si{qr), the most easily accessible quantity is 
the “equatorial” cut through reciprocal space ^eq = 
S{qr, = 0) = S'o+2S'i obtained under normal incidence 
(a = 0 in Figured]). We therefore subtract ^eq from the 
re-gridded data (an example is shown in Figure |3l left) 
and fit it with an appropriately modified version of Equa¬ 
tion O- The best fit for the data in Figure [3] is given in 
the right image. 

The vertical streak on the right in the experimental 
data is due to an imperfectly subtracted diffuse scatter- 



FIG. 4. Partial structure factors for the sample with n = 
1.97610“® A extracted from the data in Figured) Sym¬ 
bols: the equatorial structure factor Seq (experimental data). 
Lines: the structure factors So and Si. For the solid line 
represents the estimate extracted from the configuration with 
a = 60 °, while the shaded area around the curve represents 
the uncertainty (see text). 


ing ring due to a kapton window in the beam (we found 
that proper background subtraction is quite difficult for 
tilted samples, i.e. for a ^ 0). Otherwise, the agree¬ 
ment between data and the fit is good, showing that the 
model is accurate (in particular, there is no need to in¬ 
clude higher-order partial structure factors). 
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This treatment yields Si{qr), which is then subtracted 
from the equatorial structure factor Seq{qr) to yield 
So{qr). Figure 0] shows these functions for the most con¬ 
centrated sample (extracted from the data in Figure 0]). 
In the case of Si, the solid line represents the best fit to 
scattering data taken for a = 60° (Figure (21), while the 
shaded area contains the curves obtained for a = 45, 60 
and 75 providing a rough estimate of the experimental 
uncertainty. 


C. Interaction within the layer 


Once the experimental data for the partial structure 
factors So{qr) and Si{qz) (or, equivalently, Seq and Si) 
is available, it can be described in terms of the interaction 
potentials Vo(r) and Vi (r) via well-known methods in the 
theory of liquids, presented in Appendix [BJ The main 
goal of this work is determining the membrane-mediated 
potential Vo{r). The interbilayer interaction Vi{r), al¬ 
though needed for a complete analysis, is inessential and, 
furthermore, is probably very sensitive to the swelling of 
the phase. We therefore use a single functional form for 
it, given by a linear decrease from a maximum Ui when 
the particles have the same in-plane position r to 0 when 
they are a distance away. After extensive tests we 
found that, for various shapes of the in-plane potential 
Vo(r), the best fit for Vi(r) is described by an amplitude 
(7i around 1 and a range of about 25 A. 

In the following, we concentrate on the parameters de¬ 
scribing Vo(r). We considered several functional forms, 
such as an exponential or a linear decrease, a Gaussian 
centered at the origin or at the contact and a comple¬ 
mentary error function. The best result for each of them 
is shown in Figure [SJ except for the linear decrease which 
yields much worse fits. 

We obtained very good results for the Gaussian cen¬ 


tered at contact: Vo(r) = Uc exp 


r-2R 


(solid 


line in Figure [5] the corresponding structure factors 
are plotted in Figure 121 along with the experimental 
data) with fit parameters Uc = 2.2 ± 0.2 keT and ^ = 
10.8±0.8 A while the hard disk radius R is fixed at 4.7 A. 
The shaded area corresponds to the fit uncertainty, esti¬ 
mated as in Appendix IB 31 


D. Hydrophobic matching 

In the treatment above, we considered that the inter¬ 
action potential is independent of the particle concentra¬ 
tion. However, one can infer from a simple elastic theory 
(and confirm by experimental investigations^^) that the 
amplitude of the interaction decreases with the concen¬ 
tration of inclusions; briefly, inserting a new particle is 
easier when the bilayer is already deformed by existing 
inclusions. Accounting for this effect requires introduc- 


3.0 -1 

2.5 - 

2.0 - 

H 

m 

1.5 - 

1.0 - 

0.5 - 


0 



FIG. 5. Best results for some of the functional forms used for 
Vo{r) and corresponding values. For the Gaussian centered 
at contact, we also plot the uncertainty range obtained as 
described in Appendix IB 31 


ing a new fit parameter, the concentration no, for which 
the interaction decreases significantly. We will consider 
the simplest model, whereby the membrane deformation 
is linear (and thus the elastic energy is quadratic) in the 
concentration: Vb(n,r) = Vb(0,r)(l — n/ng)^ for n < ng 
and Vb(n,r) = 0 for n > ng. We take Vi(r) as inde¬ 
pendent of the concentration. As before, the in-plane 
potential is described by a Gaussian centered at contact. 
The fit is significantly improved, yielding = 1.00, for 
parameters Uc = 5.0 ± 0.6 keT, ^ = 9.3 ± 0.7 A and 
no — 5 10“^A (roughly twice the highest experimen¬ 
tal value available). 

We emphasize that the values presented above corre¬ 
spond to the highly diluted limit Vg(0, r), i.e. to the inter¬ 
action between two particles in the absence of any other; 
this is the relevant case for comparison with most theo¬ 
retical and numerical models. 

One can see this decrease of the interaction potential 
between two particles due to the concentration (i.e. to 
the presence of other particles) as a crude way of ac¬ 
counting for many-body effects. However, this approx¬ 
imation is not controlled and its accuracy can only be 
verified by more complete theoretical models or by nu¬ 
merical simulations^. 


V. DISCUSSION AND CONCLUSION 

The main results of this work are the in-plane interac¬ 
tion potential of the BuSnl2 particles, Vg(r) in the lamel¬ 
lar phase and V (r) in solution; they are summarized in 
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FIG. 6. Best fits (lines) to the experimental data (symbols) with the model described in Section [IV Cl Vo(r) is a Gaussian 
centered at contact (solid line in Figure [5]) 


Figure [T] One can see straightaway that the uncertainty 
is small compared to the amplitude of the curves and to 
the difference between them, except very close to contact. 

The qualitative conclusions to be drawn are similar to 
those already obtained in Reference [13 using a simpli¬ 
fied model, although the amplitudes are somewhat dif¬ 
ferent. The interaction potential in solution is short- 
ranged (about 5 A from contact), but at short distances 
it becomes stronger than the in-plane repulsion estimated 
with the constant-amplitude model (dfe^T, as opposed 
to 2.2 ksT at contact). Notwithstanding the uncertainty 
and the non-additivity of the potentials, it is tempting 
to speculate on a possible non-monotonic dependence of 
the “effective” membrane-mediated interaction, defined 
by subtracting the interaction in solution from that in 
the lamellar phase. 

Appendix A: Reproducibility 

Since the signals to be measured are fairly weak and 
unstructured, the question of their reproducibility is of 


the utmost importance. We therefore measured a few 
samples twice (about one year apart), in the same exper¬ 
imental conditions (see Section HD, and in normal inci¬ 
dence, so that only the equatorial structure factor Seq{qr) 
is available. The results are shown in Figure H] (left). 
Gray lines: earlier data, used in the analysis of Ref. 
Black lines: more recent data, measured on the same 
samples. For clarity, the curves are shifted vertically in 
steps of 0.25. 

Notwithstanding the discrepancy close to the peak po¬ 
sition, the agreement is quite good, especially at small 
angles where S departs significantly from 1. Indeed, this 
is the range where the interaction potential plays a signif¬ 
icant role. For coherence, the comparison is done using 
the same simplified treatment as in Ref. [U, i. e. neglect¬ 
ing the interaction between layers Vi(r) and obtaining 
the Fourier transform V{qr) of the remaining interaction 
potential Vb(r) = V{r) via the RPA approximation; see 
Ref. for the details. 

The two sets of data for V{qr) are shown in Figure [5] 
(right). The solid line and the shaded area correspond 
to the results of Ref. [13 (average value and uncertainty) 


























































































FIG. 7. Best results for the in-plane interaction potential 
Vo{r) in the lamellar phase, with constant amplitude (solid 
line - see Section IIV Cl) and with variable amplitude (dotted 
line - see Section FlV Dl) as well as for the interaction potential 
V{r) in solution (dashed line - see Section FlV Al l The shaded 
areas represent the uncertainty, determined as explained in 
Appendix IB 31 

while the symbols are obtained directly from the more 
recent values of Seq{qr) (black lines in the left panel). 
Although the new data seems to yield a stronger interac¬ 
tion than the old, the difference is within the estimated 
uncertainty range. We conclude that the measurement 
method is reliable and that the samples did not age sig¬ 
nificantly over a period of one year. 

Appendix B: Interaction model 
1. Analytical development 




FIG. 8. Left: comparison between the equatorial structure 
factor Seciiqr) determined in two different experimental runs 
(black and gray curves). A noticeable difference appears for 
n = 0.584, signaled by the double-headed arrow. Right: 
comparison between the corresponding potentials V{qr), esti¬ 
mated via the RPA approximation. 

This system has been studied for a long time^i^^. In 
the following, we will use the notations of Ref. except 
for the interaction potential, denoted here by V (instead 
of (j)) for coherence with previous work. We consider n 
bilayers (indexed by lowercase Greek indices) with peri¬ 
odic boundary conditions (bilayer n is identical to bilayer 
0). The relevant equations are: 

The definition of the total correlation function h(r) in 
terms of the radial distribution function g{r)'. 

hap{r) = gaii{r) - I (Bl) 

The Ornstein Zernike relation, defining the direct cor¬ 
relation function c{r): 

^a/3(^) — CcK^(r) -1- ^ ^ I d S Cc^..y (I s I )/i^^ ( [t* s|) 
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(B2) 

where pj = p is the number density in the plane of 
the bilayers and 7 runs over all bilayers. Alternatively, 
the equation above can be rewritten in reciprocal space, 
yielding: 


The interaction is treated in the framework of the 
Ornstein-Zernike relation, with the Percus-Yevick closure 
approximation. The structure factors Sm{q) are com¬ 
puted using the method introduced by Lado^i*^. An 
implementation for the case of no interbilayer interac¬ 
tion has already been used^^. In the present case, we 
extend the model to the case of several interacting bi¬ 
layers. As discussed in the main text, the system can 
be mapped to a multi-component fluid, where each com¬ 
ponent corresponds to particles inserted within one bi¬ 
layer and all particles formally occupy the same two- 
dimensional space. Throughout this analysis, only the 
two-body interaction is taken into account. 


hapiq) = Cap{q) + ^ p-fCa-^{q)h^p[q) 
1 


The Percus-Yevick closure: 


Cocp{r) 


[1 


Kp{r)\ 


1 — exp 


(Y^\ 

V ksT ) 


(B3) 


(B4) 


with Vap{r) representing the interaction potential be¬ 
tween particles situated in bilayers a and /3. For sim¬ 
plicity’s sake, we consider that the particles only interact 
within the same bilayer and with particles in the nearest- 
neighbor bilayers {Vapir) = 0 for |a — /3| > 1). Further¬ 
more, translation symmetry (an identical environment 
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for all bilayers, i.e. periodic boundary conditions) and 
mirror symmetry {Vap = Vsq, and similar relations for 
functions c and h) allow us to define all functions with 
respect to layer 0, i.e. a = 0 throughout. 

The equations above must be solved numerically; for 
convergence reasons, it is much better to work with 
the indirect correlation function (denoted in Ref. [3l| by 
H{r)), 70 / 3 ( 7 ’) = hop{r) — cop{r) instead of the total cor¬ 
relation function. 

As an example, for n = 3 bilayers (with periodic 
boundary conditions, 3 = 0), equations (IB3I1 and (IB4I1 
can be written out explicitly as: 

700 (g) = P [coo(7oo + Coo) + 2 coi(7oi + cqi)] (B5a) 

701 (g) = P [coo(7oi + coi) + coi(7oo + coo) + coi(7oi + coi)] 

(B5b) 

coo(7’) = [1 -b7oo][exp(-V()o/feB'r) - 1] (B5c) 

coi(r) = [1 -|-7oi][exp(-V()i/A:Br) - 1] (B5d) 

where we omitted the argument of the functions on the 
right-hand side. One starts by applying in real space 
equations (IB5b -dl with reasonable initial guesses, and 
then solving ()B5h -bl in reciprocal space. The procedure 
is iterated until stability is reached. The experimentally 
relevant structure factors are simply: 5'o(g) = 1-1- /loo(g) 
and Si{q) = hoi{q). For completeness, the solution of 
(IBfih -bl is: 

( \ _ “Coo(coo + Coi - 1/p) + 2c§i(coo + 1/p) 

~ (coo - 1/p)^ + COI (coo - 1/p) - 

(B6a) 

, s _ -(cqo + Coi)(coo - 1/p) + Cqo/p -f 

7oi g coi coi(coo -1/p)- 

(B6b) 

We checked that, for relevant values of the interaction 
parameters, the structure factors obtained with n = 5 
and n = 7 are almost superposed, and slightly different 
from those obtained for n = 3 (with periodic boundary 
conditions in all cases). All results presented above are 
obtained with n = 7. We also checked that the higher or¬ 
der structure factors are negligible, viz. S 3 {q) S 2 {q) 

S'i(q), in agreement with the experimental observations. 

2. Software implementation 

The routine discussed above is implemented in iGOR 
Pro 6 as an ” all-at-once” function, i.e. all points of the 
output vector are returned at the same time. This ap¬ 
proach is particularly useful for such iterative procedures, 
where, in order to get its value at one point, the en¬ 
tire function needs to be calculated anyway. The Fourier 
transform was implemented as a matrix operation with 
pre-coinputed coefficients over a fixed equidistant grid. 


Wherever possible, the wave operations were done using 
the MatrixOp command. 

For each combination of fit parameters, the structure 
factors So and Si are calculated for all densities. Com¬ 
parison with the experimental data yields the goodness- 
of-fit function which is minimized using the Optimize 
operation with the simulated annealing method. 

3. Confidence range 

Once a minimum is found, the goodness-of-fit function 
is plotted as a function of the parameters Uq (or Uc) 
and f. 

For the 3D case (interaction in solution), one notices 
that the minimum is in fact an extended valley, covering 
a wide range of parameters, and roughly defined as the 
locus of the points where the integral of the potential is 
constant. This is understandable, since the fit is mainly 
sensitive to the low values of the scattering wave vector 
q, where it is affected by F(g —>■ 0 ), i.e. by the integral 
of V{r), rather than by its finer details. We therefore 
choose a number of points along this valley, with close 
to minimum y^, and trace the corresponding V{r) in¬ 
stances. The confidence range (shown as a shaded area 
in FigurelU bottom) is chosen by manually adjusting the 
parameters to yield lower and higher envelopes to this 
sheaf of curves. 

In the 2D case (interaction within the layer), the min¬ 
imum is well localized in the (C4iC) plane. However, a 
standard statistical estimation of the confidence range is 
irrelevant, since it would yield extremely tight confidence 
ranges. Indeed, the discrepancy between model and data 
is mainly due to systematic errors and (presumably) to 
an inexact functional form for the interaction potential. 
In order to account for the latter effect we took as an ac¬ 
ceptable increase Ay^ the difference between the Gaus¬ 
sian centered at contact (y^ = 1.34) and that centered 
at the origin (y^ = 1.40), both plotted in Figure [51 The 
uncertainty on Uc and f given in Section II V Cl is deter¬ 
mined based on this numerical value. For comparison, we 
also show the exponential decay with y^ = 1.52. We fol¬ 
low a similar procedure for the variable amplitude model 
discussed in Section IIV PI with Ay^ = 0.1. 
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